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I. INTRODUCTION 



It often occurs that a physical problem can be treated 
by means of rather different mathematical approaches. 
Problem of pair creation from vacuum by electric field 
was first considered on the basis of the first quantized 
relativistic theory in the context of the Klein paradox, 
i.e. scattering at electrical potential barrier [lrfcj . Later 
on, this problem was studied in the framework of quan- 
tum field theory by the effective action technique 
by explicit construction of ^-matrix in the non-stationary 
gauge 04l3l and by semiclassical methods [HI [l2| , see 
also the monographes fl3l - il5j for review of these and 
some other approaches. 

More recently, there appeared yet another approach 
to pair creation effect in external field |16l420| . which is 
based on equations that bear strong resemblance to ki- 
netic equations widely used in plasma physics and other 
non-equilibrium problems. Following the authors of these 
papers, let us call it the quantum kinetic equation (QKE) 
approach. The obvious benefits of the QKEs are simple 
implementation of numeric calculations [lj], Oil H3 and 
ability of natural inclusion of account for back reaction 
to the external field 



21| . On the other hand, there are 



also several disadvantages, including (at least, currently) 
restriction to homogeneous fields only and seemingly non- 
gauge invariant form, non-ability of finding exact analyt- 
ical solutions even for the problems that knowingly ad- 
mit such solutions, and, as we are demonstrating in this 
paper, delusive physical interpretation of the basic quan- 
tities involved in QKE. Obviously, synthesis of different 
approaches may be very helpful for further development 
of the theory. 
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The QKEs had been rigourously derived from the first 
principles of QED for a specific case of homogeneous 
time-variable electric field for creation of both scalar 
[HI, H(| and fermion [T3, EH pairs. Therefore, this ap- 
proach must be manifestly equivalent to other exact 
methods mentioned above. Surprisingly, to the best of 
our knowledge, no attempts were made in the literature 
to establish this equivalence explicitly until the recent 
paper 22]. Moreover, let us mention the attempts to 
contrast the results obtained by QKE to those obtained 
by more traditional approaches 23] with finding differ- 
ence among them, and to derive the different source term 
for QKE starting from the S- matrix approach [24|. As 
for the paper [22j , though it have stated for the first time 
the aforementioned equivalence, explicit correspondence 
was provided only for scalar pair production and, in ad- 
dition, in the form which is not best suited for practical 
purposes (e.g., for testing the numerical routines with the 
exact solutions). 

In order to simplify the QKE, several ad hoc approxi- 
mations had been suggested, such as, e.g., the Markovian 
and the low-density approximation [23l . |25| . However, ex- 
pectations that these approximations can be justified un- 
der some conditions are based on the analogies with the 
properties of more usual kinetic equations. An important 
point here is that the function satisfying the QKE does 
not actually possess the meaning of particle distribution 
as long as particle creation process goes on. Rather, this 
function accounts for a "quantum soup" of real particles 
and vacuum fluctuations. These ingredients can not be 
separated in principle owing to the uncertainty relations. 

In the present paper, after a brief review of the QKE 
approach in the section|TTl we provide the explicit bilinear 
ansatz that converts QKE to the (dimensionally reduced) 
Klein-Gordon or Dirac f Sec. equation. Our ansatz al- 
lows to adopt the known exact analytical solutions of the 
pair creation problem in the QKE context. In particu- 



2 



lar, we illustrate our correspondence in more details for 
the paradigmatic case of constant electric field. In the 
Sec. IIVI we explicitly demonstrate that the low-density 
approximation is not asymptotically exact in the weak 
field limit and estimate its accuracy. The discussion and 
concluding remarks are collected in the Sec. [V] 



II. QKE APPROACH: ILLUSTRATION WITH 
ID OSCILLATOR MODEL 



Derivation of QKEs in the context of pair creation 
problem is thoroughly discussed in the literature [13, EH, 
20. 22]. However, for convenience of the reader and in 
order to introduce the notation, let us provide here a 
sketch of derivation of QKE in a more simple toy model 
problem - for parametric excitation of a ID quantum os- 
cillator with time varying frequency ui(t). The familiar 
Hamiltonian of the oscillator reads 



H 



p2 

2 



\t)Q 2 



where Q and P are the coordinate and the momentum 
operators, obeying the ordinary commutation relations. 

Consider the approximate WKB solution for the oscil- 
lator eauation|29j Q + w 2 (t)Q = 0, 



*(t) 



e- <e(t) , 9(t)= u(t)dt, (1) 



and the time-dependent operators a(t), a'{t), defined in 
such a way that the relations 

Q{t) = §{i)a{i) + $*{t)(J{t), 
P{t) = -iw{t)m)a{t)-^*{t)a\t)l (2) 

after all are carried out exactly. These operators are ob- 
viously expressed in terms of the original operators Q 
and P via 

a(t) = (u(t)Q + iP)$*(t), a t(t) = (w(t)Q - »P)#(t), 

(3) 

are time- dependent and obey the commutation relation 
[a,at] = l. 

Next, by taking into account the equations of motion 
for the oscillator, we obtain: 



a{t) 



*(*) 
2u{t) 



; 2i0 «at(t), dt(t) - ^L e -^W 



2a; (i) 



a(t). 



Thus, time dependence of the introduced operators is 
caused solely by variation of the oscillator frequency. 

Now, consider the "instant excitation number" N(t) — 
a^(t)a(t). By taking into account the Eqs. (@|, we can 
write 



N(t) 



*(*) 
2u(t) 



(V + V^), 



(5) 



where we have denoted V(t) 



-2i0(t) n 2 



a 2 (t). On the 



same grounds, the time derivative of V(t) can be ex- 
pressed in the form 



V(t) = -2iuj{t)A + 



2ui{t) 



(l + 2N(t)). 



(6) 



The averages Af(t), V(t) of the operators N(t) and V(t) 
over the actual quantum state of the oscillator obviously 
obey the same equations ((5]) and ((6]). 

Suppose that initially (i.e., at t — > — oo), the frequency 
was constant and the oscillator was hosted in some sta- 
tionary state \n). This means, in particular, that the 
operators a and had represented just the usual low- 
ering and raising operators acting on the ladder of the 
stationary states. In particular, the "anomalous average" 
V(— oo) is zero initially under our assumption. With such 
initial condition, the equation ^ can be integrated out 
with the result 
t 

V(t)= / di'i^l[l + 2AA(0]e^ [e(t) ' e( *' )1 . (7) 
J 2w(t ) 

— oo 

The analogue of the QKE for our toy model 



M{t) 



— oo 

(8) 

can be now obtained by substitution of the Eq. (0 into 
the Eq. ©. 

Let us demonstrate how the integro-diffcrcntial Eq. , 
or, equivalently, the system of differential equations (|5|). 
(j6]) can be converted back into the original oscillator 
equation. Let ^(t) be an exact, positive frequency so- 
lution of the equation 



* + LJ 2 (t)^ = 0. 



(9) 



Under our assumptions, the latter means that (up to a 
constant phase factor) ^(t) — > $(£) as t — > -co. 

The position and the momentum operators of the os- 
cillator can be expressed in terms of the aforementioned 
solution by 

Q{t) = ^(t)a+t>*(t)a\ P{t) = i>{t)a+i>*{t)a\ (10) 

where the operators a, are now time-independent. 
These operators obey the usual commutation relations 
[a, ar] — 1, provided that the solution ^(t) is normalized 
by 



(11) 



and they acquire the meaning of in-operators (i.e., de- 
fine the ladder of the initial states). By averaging the 
operator equalities 



V(t) 



2uj(t) 

! (t)Q 2 - P 2 - iuj(t)(QP + PQ) 
2uj{t) 
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which follow from the Eqs. (|3]) over the initial state |n) = 
(a t )"|0)/\/^! (where a|0) = 0), we find that 

^) = n + " a ^'^'' + 't (t)|a - a;(f) (l + 2n), (12) 



2w(t) 



and 



V(*) - I - 2*I*WI ) + 2n >" 

(13) 

It is simple to prove by direct substitution that the formu- 
las (fT2|) . ([T3|) indeed define the exact solution of equations 
([5]), © with initial conditions Af(— oo) = n, V(— oo) = 0, 
and hence, of equation (|SJ), see also the Appendix 1X1 



where is defined by the same expression (|T5|) . 

The QKE for the fermion case differs in its form from 
the Eq. © [rEEl, 

— oo 

xcos[2(e(t)-e(t'))] ) (is) 

at the same time the analogues for ([5]), ((6]) read 



(19) 



III. SCALAR AND FERMION PAIR 
PRODUCTION 

Consider a charged scalar field ip affected by a linearly 
polarized time-dependent homogeneous electric field, di- 
rected along the z axis. In this section, we assume for 
simplicity that the quantum field is initially in its vacuum 
state. After separation of variables, the Klein-Gordon 
equation 

[£>„£>" + m 2 ] <p(t, r) =0, £>„ = d M + icA^ (14) 

for the function tp(t,r) = (2-K)~ 3 ^ 2 ^ p (t)e lpT is reduced 
to the form with 

u(t) = \Ar 2 + e 2 , n = Pz + eA(t), e 2 =p 2 ± + m 2 . 

Although this problem deals in fact with a complex, 
rather than the real valued oscillator, nevertheless the 
QKE acquires the same form ([5]) as above [l8|, |2(J. Thus, 
its solution must be also of the same form (fT2| , (fl3)l as for 
the real valued oscillator. In accordance with our initial 
conditions, we only have n p = 0. 

The actual difference between the two cases is only in 
the meaning of the quantities Af and V. Namely, if a p and 
b p are the destruction operators for particles (with charge 
— e) and antiparticles (with charge e), respectively, then 
Af p = (vac | a p a p | vac) is the "instant average number of 
particles" (or pairs), whereas the "anomalous average" 
is now defined as V p = e~ 2iep O(vac|a p 6p|vac). Note 
that these "anomalous averages" are also important in 
condensed matter problems, e.g. they participate in for- 
mulation of semiconductor Bloch equations [28| . 

Let's now turn to the case of the fermion (spinor) field. 
The Dirac equation 

[i 7 "£>„-rr#(i,r) = 0, (16) 

after transition to its squared version [D^D^ + m 2 — 
i(e/2)F tlu j fl ^]tp = 0, for the function of the form 

ip(t,r) = (2n)- 3 / 2 ^ ( p ± \t)e ipr u ± (7°7 3 u± = ±u±) reads 
+ [w 2 (t)±ieA(t)]y( ± ) = 0, (17) 



V(t) = -2iw(t)V + ^-±±(1 - 2Af). (20) 

The functions Af and V are defined similarly as for the 
scalar field, but are carrying an additional two-valued 
spin index. We omit this index as well as the momentum 
index for simplicity, note that the coefficients of the QKE 
are independent on the spin index. 

In the fermion case, reduction of the equations (|19l) . 
(f2"0"]) to the Eq. (JTTJ) is less trivial than in scalar case, 
mostly due to a different "statistical weight" factor 1 — 
2Af on the RHS. However, this procedure can still be 
done, e.g. following the derivation of the QKE. Let us 
present here only the result, which can be easily proved 
by direct substitution. Namely, we have 

^(*) = TT-A tHI* (±) | 2 + l* (±) | 2 

lujyu) — 7TJ 

Ti^ (±) ^> {±) * -*( ± )*( ± )*)}, (21) 

v(t) = -— fw 2 !*^! 2 -!*^! 2 

^(* (±) ^ (±) * + ^ (±) * (±) *)}, (22) 

where the functions obey the Eq. (TT71) with upper 

and lower sign, respectively, and are subject to the nor- 
malization condition 

w 2|^(±)|2 + |^(±)|2 T in{¥ ±)y(±)* _ ^(±)^(±)*) = £ 2 ; 

(23) 

see Appendix [B] for more details. In practice, either of 
the functions <3/(±) can be used to compute the quantities 
Af and V. 

In the rest of the section, let us apply our results to the 
particular case of a constant electric field, A(t) = —Et. 
As it was stressed in the introductory section, though this 
problem is known to be solved exactly in terms of the 
parabolic cylinder functions _D„, no exact solution had 
been ever presented previously for QKEs. The positive 
frequency solutions for this case are explicitly given by 
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FIG. 1: Occupation number M (left panel) and the magnitude of "anomalous average" V (right panel) for scalar (dotted line) 
and fermion (solid line) fields versus the dimensionless time — £ = s/eE(t — p z /eE) for A = 2.5. The dash lines show the 
asymptotic values at t — > +oo. 



(compare to 0413) 



^scal(t) 



3 -ttA/8 



[2eEy/ 



—D^[{l-i)Z], 



(-) 

ferm 



(t) 



(24) 



(25) 



where A = e 2 /eE and £ = —\feE(t — p z /eE). These 
formulas, together with the general expressions (IT2l . (p~3|) 
and (|2"T1) . (|2"2"j) . define the desired solution for QKE in 
constant electric field. 

The dependence of the quantities AA and |V| on time 
for both cases of scalar and fermion field is shown at 
Figs. Q] and [2] for A = 2.5 (subcritical pair production) 
and A = 0.1 (supercritical pair production), respectively. 
The initial state is always assumed to be vacuum. The 
plots JV(t) obtained by direct numerical integration of 
QKE can be found on Figs. 2 and 3 of the Ref. (2j|. It is 
clear from the figures, that the curves are starting from 
zero at t — > — oo in accordance with our initial conditions 
and are tending to some asymptotic values at t — > +oo. 
These asymptotic values are AA(+oo) = e~ 7rA (for both 
scalars and fermions) and |V(+oo)| = e~ 7rA / 2 (l±e~ 7rA ) 1 / 2 
(for bosons and fermions, respectively) and agree with 
the results obtained by other, e.g., S-matrix, methods. 

However, this transition is not monotonous but is ac- 
companied with oscillations in some transient region of 
width At ~ y/X/eE around t = p z /eE, which in fact es- 
sentially coincides to the commonly accepted "coherence 
length" (or formation time) 0413 f° r a P a i r production 
process [30]. It can be also observed from the figures that 
the magnitude of oscillations is relatively much larger 
(compared to the corresponding asymptotic values) in 
the case of large values of A, i.e. for subcritical pair 
production. These oscillations can not be attributed to 



instant increase in particle production itself, because the 
notion of particles and antiparticles can not be defined 
rigorously in the transient region. This must be clear 
already from the N — N commutation relation 



[N,N] 



(v-vi) 



(26) 



which follows from the Eq. ([5]). As a consequence, we 
have the uncertainty relation 

|wImV| 



AAA- AAA > 



(27) 



Similar uncertainty relation can be derived for fermions. 
In virtue of (|27|) . oscillations in the transient region are 
accompanied with quantum fluctuations. Similarly to 
AC, the quantity ImV on the RHS also oscillates with 
frequency 2w. The quantity AA characterizes simultane- 
ously the real particles and the vacuum perturbations 
and there is no way to extract information relevant for 
the particles alone. Hence, it is rather more accurate 
to use the term "quasiparticles" instead of "particles" in 
this context. Only outside the transient region, where 
ui varies slowly, quasiparticles turn into the real parti- 
cles. This means that, in contradiction with the claims 
of some other authors, particle creation process can not 
be resolved on the temporal and spatial scales of the co- 
herence length. We'll come back and give more comments 
on this peculiarity in the sequel. 

The total number of pairs created in a unit volume 
during the time period t\ < t < t% can be obtained as 
usually by summing over all the excited modes of the 
field. For scalar field, this quantity can be computed as 



d 3 p 



■El-, 



A/p(+oo) = 



dp z f d 2 p± 



eEti 



(2tt) J (2tt) 2 



e 2 £ 2 
8tt 3 



■(*2-ti), (28) 



5 



1,0-1 1 ' 1 ' 1 1 1 1 1 1 1 ' 1 1 1 I 1,4i ' ' ' ' ' ' ' ' ' ' ' ' ' ' L 




FIG. 2: Occupation number M (left panel) and the magnitude of "anomalous average" V (right panel) for scalar (dotted line) 
and fermion (solid line) fields versus the dimensionless time — £ = yeE(t — p z /eE) for A = 0.1. The dash lines show the 
asymptotic values at t — > +oo. 



where we have taken into account that only particles with 
ti < Pz/ e E < t 2 are being created in the time interval un- 
der consideration. For fermion case, the pair production 
yield is twice larger due to additional summation over the 
spin index of the modes. These results are well known, of 
course, and had been obtained previously by other meth- 
ods, though in some modern papers the pair production 
yield is sometimes incorrectly identified with the doubled 
imaginary part of the Heisenberg-Euler-Schwinger effec- 
tive action. The latter is related in fact to the vacuum- 
vacuum transition probability and has no exact relation 
to the pair production yield. 



approximation does not provide good accuracy (though, 
of course, gives the correct order) in the whole range 
of the field strength, whereas it is believed that the low- 
density approximation works good for weak fields. In this 
section, we evaluate the total pair production rate within 
the low-density approximation analytically for weak con- 
stant field and explicitly estimate its accuracy. 

For definiteness, let us illustrate the line of computa- 
tions using the Eq. ([5} and later come back to the fermion 
case. Within the low-density approximation, the total 
pair production rate per unit volume is defined by 



IV. LOW-DENSITY APPROXIMATION 

The QKEs ©, (|T8)l look very complicated, in partic- 
ular they possess non-Marcovian character. The latter 
means that the RHS of the QKE depends on the history. 
Of course, this dependence arises just for elimination of 
the equations flBJ), (12TJ1) . so that the problem admits an 
equivalent formulation in local terms. Nevertheless, it 
may be useful in some problems to deal with QKE on 
its own footing, without reference to differential equa- 
tions that lie in its origin. Several attempts were made 
in order to simplify the QKE, including the Markovian 
(taking the factor 1 ± 2AT outside the sign of the integral 
over t') and the low-density (neglecting Af on the RHS) 
approximations [H, [25| . Intuitively, they would be rea- 
sonable if AT was large and monotonously increasing, or 
essentially small, respectively. 

However, the occupation numbers are actually not 
large and, besides, the oscillatory behavior discussed in 
the preceding section generally contradicts both assump- 
tions. In practice, it was shown by direct numerical in- 
tegration in the papers cited above, that the Markovian 



R ld = 



d 3 p U! p (t) 

(2tt) 3 2uj p (t) 



dt' ■ 



>(*') 



>(*') 



cos 



u p (t") dt' 1 



(29) 

where we have recovered the momentum index of the 
modes. For the case of constant field, we have ir(t) = 
p z — eEt, so that we can change the variable of inte- 
gration p z by w — ir(t)/e. The time derivative io p {t) = 
eEir(t)/u} p {t). Also, let us come to the new variables 
u = n(t')/e = [n(t) + eE(t - t')]/e and v = n(t")/e = 
[w(t) + eE(t — t")]/e instead of t' and t", respectively. 
After these transformations, we are coming to 



eE f d 2 p A 



dw 1 



du i 



(27r) a J 1 + w 2 J 1 

— CO w 

x cos {2\[g{u) - g(w))} , 



(30) 



where g{x) = f Q VT + v 2 dv. The integrand is invariant 
under exchange ufjiu, therefore it is possible to repre- 
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sent the expression on the RHS in the form 

oo 

— OO 

The same representation is valid in the fermion case 
as well, but with an additional factor two due to sum- 
mation over the spin variables and with ff e rm( x ) — 
f^due^W/il + u 2 ). 

In the weak field limit A is large, so that the functions 
fscai, fferm can be evaluated by the approximate meth- 
ods (see the Appendix [C]). Thus we obtain 

/ SCQ K2A)-^ e - A / 2 [l + 0(A- 2 / 3 )], 

/ /erm (2A) = ye-" A/2 [l + 0(A- 2 / 3 )]. (32) 

The integral in (13"TT) is evaluated easily and we discover 
that even for A 3> 1 the total pair production rate in 
the low-density approximation exceeds the exact value 
given by (gHJ by the factor (tt/3) 2 , i.e. by 9.6%. We 
attribute this discrepancy to contribution of oscillations 
in the course of pair creation that have been discussed 
previously. 

V. DISCUSSION 



computations. This is not the case, however, since the 
QKE approach admits simple incorporation of backreac- 
tion on the external field [21(. In addition, the form of 
the kinetic equation is very pleasant in principle for in- 
corporation of the secondary processes in the resulting 
electron-positron plasma [26l . |27[, such as annihilation, 
hard photon emission, Compton scattering, etc. by sim- 
ple inclusion of the appropriate collision terms on the 
RHS. 

At the same time, in our opinion, at the present stage 
of development of QKE approach the latter attempts are 
too hasty, at least in the high-intensity laser problems. 
The first reason is that the spatial and temporal scales 
of variation of the laser field are of the same order, so 
that spatial variation of the field should be taken into 
account as well. Up to now, this has never been done 
starting from the first principles of QED, though some 
phenomenological conjectures had been given. The sec- 
ond reason is delusive physical interpretation of the oc- 
cupation number participating QKEs at the time scales 
of the order of the formation time of pair creation pro- 
cess. As we have shown in the paper, the situation here is 
quite similar to the attempts to resolve the trajectory of a 
particle with the accuracy exceeding the limits imposed 
by the uncertainty principle. We believe that the true 
kinetic equations for pair creation, which can be treated 
by the commonly accepted methods, must contain some 
sort of averaging over the fine scale of formation for the 
pair creation process. 



Wc have demonstrated that the quantum kinetic equa- 
tion (QKE) approach to pair production problem is 
equivalent to the known for a long time S'-matrix ap- 
proach by establishing explicit correspondence between 
them for both scalar and fermion fields. This equiva- 
lence has been already declared recently [22| , though not 
proved for the fermion case. As an illustration, we have 
applied our correspondence to obtain the exact solution 
for both scalar and fermion QKEs in the constant field. 
Of course, exact analytical solutions for QKEs can be 
now constructed for other situations that are known to 
be exactly solvable by conventional methods, e.g. for 
E(t) = E /cosh 2 {Kt). 

At the same time, any attempt to perturb the KQEs 
by applying some kind of approximation procedures (e.g., 
by imposing the Markovian, or the low-density approx- 
imations) inevitably destruct the correspondence. We 
have demonstrated it on the example of the low-density 
approximation which was believed to be well justified in 
the weak field limit. In contradiction with these expecta- 
tions, we have proved that for the case of constant field 
its accuracy all the same tends to 10% even under the 
optimal conditions. 

It could seem that, being equivalent to the conven- 
tional approaches and looking more complicated, the 
QKEs are useless for pair creation problems with the pos- 
sible exception of arranging the algorithms for numerical 
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Appendix A: Direct proof of the ansatz 

Let us show the way how our solution (fT2"j). (fT3"|) can 
be deduced directly from the Eqs. (J5j , d©]) , without any 
reference to derivation of the QKE. For this purpose, let 
us introduce the notation V = 1Z + ij for the real and 
imaginary parts of the "anomalous average" V. In this 
notation, the equations ([5]) and ((6|) take the form 

N=-K, (Al) 

K = -2uJ + ^ (V+ , (A2) 

j = 2u1l. (A3) 

Let us consider a sum of Eqs. (|A1|) . (|A2I) and introduce 
the new function T defined by N +TZ + 1/2 = Tlu. Then 
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we immediately get 

J = -T/2 
and after that in virtue of Eq. (|A3[) 

K = -f /4w, 

so that 

A/" = wT + T/4w- 1/2. 



(A4) 



(A5) 



(A6) 



Thus, all the functions participating the system of equa- 
tions (jAll) . (|A2|) . (|A3|) are now expressed in terms of a 
single unknown function T. 

Combining all the results (|A4I) , (| A5|) and (IA6|) in either 
Eq. (|A1[) or Eq. (IA2[) . we are coming to the third order 
equation for T: T'" = — 4w(wT)'. Multiplying it by T 
and integrating by taking account of the initial conditions 
T(-oo) = (ra + l/2)/w(-oo), t(-oo) = f (-00) = and 
the identity TT'" = (TT")' — T'T", we are arriving to 



T 2 

TT - — 



2lo 2 T 2 = 



(l + 2n) 2 



Now, substituting 



T = |*| 2 (1 + 2n) 



(A7) 



after simple transformations we are finally coming to 
|*| 2 |(f + cj 2 *)** + (** + w 2 **)*| 

(A8) 



1 + (\j r \E r * — vp*ij/) 2 



Since the derivative of the RHS of Eq. 
resented in the form 



can be rep- 



_ x [(# + w 2 *)** - (#* + w 2 **)*], 

it is clear that the sought-for function ^ satisfies the 
oscillator equation © together with the normalization 
condition[31J pip , i.e. is a positive frequency solution at 
i->-oo. Finally, substituting ([A"7]) into the Eqs. (|A4]| . 
(|A5|) and (IA6|) and using ([9|) to eliminate the second- 
order derivatives, we obtain the desired expressions (|12p 
and (fn 



Appendix B: Normalization condition for fermion 
problem 



In order to derive the condition (1231) . let us first mul- 
tiply the Eq. (J!?) and its conjugate by vj/(±)* and tW, 
respectively, and take a sum of both of them. The result- 
ing equation can be written in the form 



(u;2|*(±)|2 + |^(±)|2) =2UJU\¥ ± '>\ 2 

±in(¥ ± ^¥ ± > - ^(±)<jr(±)*). (Bl) 



Next, multiplying the Eq. (flT)) and its conjugate by 
xJ/W* and respectively, and taking their difference, 

we obtain 

(#(±)*(±)* - ijr(±)#(±)*)' ± 2i 7 r|*( ± )| 2 = 0. (B2) 

Taking into account that Co = tttt/lj and expressing 
2 7r |4i( ± )| 2 from (|B2|) . we can integrate the Eq. (|BT1) with 
the result 

w 2|^(±)|2 + |^(±)|2 Ti7r ^(±) $ (±)._^(±)^(±)* ) = congt 

(B3) 

The constant on the RHS can be evaluated in the region 
where uj(t) varies slowly so that WKB approximation can 
be applied to the Eq. (TlTl) . For this purpose, we restore 
the Planck constant in this equation as follows, 

fr^W + [u; 2 (i) ± i/wr(t)]tf (±) = 0, (B4) 

and seek for the solution in the form VfrW = e lS / h , where 
S = Sq + hSi + 0(h 2 ). Substituting this ansatz into 
Eq. (IB4[) and equating the coefficients at the same pow- 
ers of h, after all we find the WKB positive frequency 
solution in the form 



(±) 



C± 



("W±^) 1/2 e-jf«Mt. (B5 ) 



Being the components of a unique spinor, vp(±) are nor- 
malized by |*(+)| 2 + |*(-)| 2 = 1, so that C± = l/y/2. 
Thus, assuming K=l, substituting (|B5I) into the LHS of 
the condition (IB3|) and taking into account that in the 
framework of WKB approximation 5'( ± ' ^2 — icj^^ we 
are coming to the Eq. 



Appendix C: Calculation of the functions f sca i and 

fferm 

Both functions are calculated in a similar manner. 
Let us start by evaluating the function ff e rm( x ) = 
duF(n,x), F(u,x) — e ixg ^ /(1+tt 2 ), assuming that 
x is large. The stationary phase method can not be ap- 
plied in our case, because the integral has poles at the 
stationary points u = ±i of the exponent and, in addi- 
tion, these points are the branch points for g'(u) rather 
than simple zeros. 

Since the function g(u) is odd, we can write 



fferm{x) = 2 Re du 



1 



Next, we have Im g(w) — > +00 as |u| — > 00 in the sector 
< argu < 7r/2. Thus we can safely rotate the ray of 
integration by the angle 7r/2, superposing it with imag- 
inary axis. In this way, the integration contour passes 
by the aforementioned point u — i from the right side. 
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Accordingly, we have 



fferm(x) = —2 Im in Res M =j F(u, X 

oo 

+ V.p. / du F(iu,x) 



If x is large as assumed, then the main contribution to the 
latter integral comes from the region of small it. Hence 
in the first approximation we can neglect the squares of u 
and v. Finally, by transition to a new integration variable 
C = Jq V^vdv = (2-\/2/3)u 3 / 2 , we are coming to 



Obviously, Res„ =i F{u,x) = — (i/2)e~ 7rx / 4 . Since the 
function F(iu,x) is real valued for < u < 1, we 
can replace the second term in the brackets simply by 
duF{iu,x). For u > 1, we can write g(iu) = 

g(i) — J" Vw 2 — 1 dv, where g(i) — in/A. Thus, after 
shifting the integration variables as u — >■ u+X, v — > v + 1, 
we obtain 



/ 



duF(iu, x) = —e 



c/4 



du 



2u + u 2 



a 



< cxp I — ix J \/2v + v 2 dv 



1 o 



Hence, we have ff e rm( x ) = (27r/3)e 7rx / 4 . In princi- 
ple, in this way the leading corrections of the order of 



0(x 



-2/3 p -ttx/4 



) can be derived as well. In the expres- 



sion for f sca i(x) the integrand differs only by the addi- 
tional factor u. However, as it follows from our deriva- 
tion, in the leading approximation the non-singular part 
of preexponent is all the same evaluated at u = i. Thus, 
f scal (x) = (27ri/3)e-^/ 4 . 
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